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Scaling  of  quantum  capacitances  and  valence  electron  detachment  energies  is  studied  for  icosa- 
hedral  and  nonicosahedral  fullerenes.  Scaling  trends  are  considered  from  zero  to  infinite  average 
radius,  where  a  fullerene’s  local  surface  properties  are  similar  to  those  of  graphene.  Detailed  density 
functional  theory  calculations  are  performed  to  determine  the  geometries  and  detachment  energies  of 
icosahedral  fullerenes,  while  values  of  these  quantities  are  obtained  for  nonicosahedral  species  from 
previously  published  experimental  results.  Strongly  linear,  quasiclassical  scaling  versus  average  radii 
rn  is  seen  for  the  quantum  capacitances,  but  on  two  different  scaling  lines  for  icosahedral  and  non¬ 
icosahedral  species,  respectively.  By  contrast,  nonclassical,  nonlinear  scaling  versus  1  /f„  is  seen  for 
the  electron  detachment  energies — i.e.,  the  valence  ionization  potentials  and  electron  affinities.  This 
nonlinearity  is  not  accounted  for  by  classical  theories  that  are  used  to  explain  trends  in  electronic 
properties  of  fullerenes  and  usually  give  accurate  quantitative  estimates.  Instead,  simple  quantum 
equations  are  derived  to  account  for  nonlinearities  in  the  metal-particle-like  electron  detachment 
energy  scaling  and  to  show  that  these  are  responsible  for  nonclassical,  nonzero  intercepts  in  the  ca¬ 
pacitance  scaling  lines  of  the  fullerenes.  Lastly,  it  is  found  that  points  representing  the  carbon  atom 
and  the  graphene  limit  lie  on  scaling  lines  for  icosahedral  fullerenes,  so  their  quantum  capacitances 
and  their  detachment  energies  scale  smoothly  from  one  C  atom,  through  Cgo,  to  graphene. 


I.  INTRODUCTION 

Fullerene  molecules  have  fascinating  and  potentially 
valuable  properties  [1].  Because  they  are  large  molecules, 
though,  containing  many  atoms  and  electrons,  detailed, 
accurate  quantum  theoretical  treatments  are  difficult  and 
computationally  intensive.  Measurements  of  the  elec¬ 
tronic  properties  of  the  isolated  molecules  are  similarly 
challenging  [2-4].  On  the  other  hand,  the  regular,  qua- 
sispherical  structures  of  these  pure  carbon  systems  [5], 
along  with  results  from  quantum  scaling  investigations  on 
other  systems  [6-9],  suggest  that  the  fullerenes’  quantum 
properties  might  vary  or  “scale”  with  their  dimensions  in 
very  regular  ways.  Understanding  this  quantum  scaling, 
as  well  as  how  it  differs  from  classically  expected  [10,  11] 
trends,  might  simplify  greatly  the  accurate  prediction 
of  quantum  mechanical  properties  for  even  very  large 
fullerenes.  Such  insights  also  might  provide  an  archetype 
for  understanding  and  predicting  trends  in  the  quantum 
electronic  properties  for  homologous  series  of  other  large 
molecules  or  nanostructures. 

To  those  ends,  here  we  explore  the  scaling  of  the 
fullerenes’  quantum  capacitances  and  valence  electron  de¬ 
tachment  energies.  The  exploration  follows  the  scaling 
trends,  starting  with  a  single  carbon  atom  Ci  and  con¬ 
tinuing  up  to  the  limit  of  infinite  average  radius,  where 
the  local  surface  structure  and  geometry  of  a  fullerene 
are  very  similar  to  that  of  graphene.  These  trends  are 
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followed  for  fullerenes  of  both  icosahedral  and  nonicosa¬ 
hedral  symmetries. 

As  seen  in  Fig.  1,  two  different,  strongly  linear,  quasi¬ 
classical  scaling  trends  are  found  for  the  icosahedral  and 
nonicosahedral  fullerenes,  respectively,  when  their  quan¬ 
tum  capacitances  [12,  13] 

Cn  =  l/(In-A.n),  (1) 

are  plotted  versus  their  average  radii  fn.  Above,  In  and 
An  are  the  first  ionization  potential  and  first  electron 
affinity,  respectively,  for  an  A-electron,  n-carbon  neutral 
fullerene  C„.  These  valence  electron  detachment  ener¬ 
gies  are  measured  in  eV,  so  values  of  quantum  capaci¬ 
tances  Cn  are  calculated  in  the  molecular-scale  units  [6,  7] 
of  fundamental  positive  charges  per  Volt  (+e/V).  Lin¬ 
ear,  symmetry-dependent  scaling  of  Cn  versus  rn  for  the 
fullerenes,  as  seen  in  Fig.  1,  is  similar  to  that  observed 
previously  for  atoms  [6,  9]  and  small  molecules  [7,  8]. 

In  Figs.  2  and  3  it  is  observed,  however,  that  the  linear 
scaling  behaviors  of  Cn  versus  rn  for  the  fullerenes  are 
produced  by  values  of  In  and  An  that  each  exhibit  dis¬ 
tinctly  nonlinear  scaling  versus  l/rn.  In  addition,  these 
figures  show  that  different  types  of  nonlinearities  are  ex¬ 
hibited  by  the  detachment  energies  of  both  the  icosa¬ 
hedral  and  nonicosahedral  species.  Figure  2  further  il¬ 
lustrates  that  such  nonlinear  behaviors  are  contrary  to 
expectations  (dotted  lines)  from  purely  classical  theo¬ 
ries  [10,  11]  that  have  been  applied  [4]  to  explain  the 
scaling  of  fullerene  detachment  energies  and  their  limit¬ 
ing  behavior  as  l/rn  approaches  zero. 

At  the  same  time,  other  nonclassical  behaviors  are  seen 
in  Fig.  1.  These  are  the  nonzero  capacitance  intercepts 
associated  with  the  quantum  capacitance  scaling  lines 
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that  otherwise  are  classical  in  form. 

Moreover,  it  is  discovered  and  proven  algebraically  in 
this  work  that  the  nonclassical,  nonzero  curvatures  in 
the  detachment  energy  trends  for  the  fullerenes  lead, 
via  Eq.  (1),  to  the  nonclassical,  nonzero  intercepts  in 
their  otherwise  classical  capacitance  scaling  lines.  These 
capacitance  intercepts  have  two  different  signs  for  the 
icosahedral  and  nonicosahedral  species,  which  correspond 
to  the  two  different  sets  of  curvatures  for  their  detach¬ 
ment  energy  trends.  This  correspondence  is  summarized 
schematically  in  Fig.  4. 

For  the  icosahedral  fullerenes,  the  aforementioned  re¬ 
sults  are  derived  from  detailed  density  functional  the¬ 
ory  (DFT)  calculations  performed  in  this  work  to  deter¬ 
mine  the  valence  electron  detachment  energies  and  aver¬ 
age  radii.  We  calculated  An,  and  rn  for  nine  such 
molecules  containing  numbers  of  C  atoms  ranging  from 
n  =  60  to  2160.  Also,  for  ten  nonicosahedral  fullerenes 
of  C2  and  T>2  symmetries,  we  tabulated  from  the  litera¬ 
ture  experimental  values  [2-4]  of  In  and  An,  as  well  as 
theoretical  values  [14,  15]  of  rn.  All  these  data,  plus  the 
associated  values  of  Cn,  appear  in  Table  I  and  are  used 
to  plot  the  graphs  in  Figs.  1,  2,  and  3. 


II.  METHOD  AND  RESULTS 

To  obtain  the  detachment  energies  for  the  nine  icosahe¬ 
dral  fullerenes,  we  performed  DFT  calculations,  using  a 
specially  developed  computer  program  [16] .  The  neutral- 
state  geometry  of  each  structure  C„  was  optimized  and 
its  average  radius  fn  determined  using  a  Gaussian  6- 
311G**  basis  and  a  functional  [17]  that  gives  an  exact 
geometry  for  Cgo-  The  DFT  calculations  involved  the 
reoptimization  of  previously  determined  [18]  icosahedral 
geometries,  as  well  as  the  optimization  of  a  new  geometry 
for  C32o-  Unlike  the  other  icosahedral  fullerenes,  which 
have  closed-shell  neutral  ground  states,  C320  has  an  open- 
shell,  triplet  electronic  state,  with  two  unpaired  electrons 
in  the  33gu  orbital.  Though  this  ionizes  to  a  doublet 
cationic  state,  as  do  the  other  icosahedral  fullerenes,  af¬ 
ter  electron  attachment  it  does  not  yield  the  usual  dou¬ 
blet,  but  forms  a  quartet  anion  with  three  electrons  in 
the  33 gu  orbital. 

After  determination  of  the  neutral  geometries,  to  eval¬ 
uate  the  electron  detachment  energies,  In=[En(N  —  1)  — 
En(N)]  and  An  =  [En(N)  —  En(N  +  1)]  for  each  n,  the 
same  program  and  basis,  but  a  different  functional  that 
gives  exact  atomic  energies  [16],  were  used  to  separately 
determine  the  total  energies,  En(N ),  En(N  —  1),  and 
En(N  +  1),  for  the  lowest  energy  neutral,  cationic  and 
anionic  states,  respectively.  Numerical  results  calculated 
in  this  way  for  the  icosahedral  fullerenes  appear  in  part 
A  of  Table  I. 

Part  B  of  Table  I  presents  the  corresponding  data 
for  nonicosahedral  fullerenes.  These  values  were  assem¬ 
bled  from  the  experimental  literature,  as  described  in 
the  previous  section.  Most  of  the  experimental  valence 
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FIG.  1.  Quantum  capacitances  Cn  of  n-carbon  fullerenes  plot¬ 
ted  versus  their  average  radii  rn.  The  solid,  dashed,  and 
dotted  regression  lines,  respectively,  in  the  main  portion  of 
the  figure,  fit  to  high  degree  of  confidence  (rn,Gn)  points 
for  icosahedral  (Xh)  fullerenes  (crosses),  for  nonicosahedral 
fullerenes  of  C2  and  T> 2  symmetries  (dark  squares),  and  for 
a  classical  sphere  model  (open  circles)  of  the  icosahedral 
fullerene  capacitors.  Regression  equations  and  parameters  for 
scaling  lines  are  displayed  in  the  graph.  The  inset  expands 
the  scale  for  the  plot  of  the  nonicosahedral  fullerene  points, 
their  scaling  line,  and  the  classical  scaling  line.  See  text. 


electron  detachment  energies  available  in  the  literature 
for  nonicosahedral  fullerenes  were  for  species  of  C2  and 
T> 2  symmetries.  Limited  data  [3,  4,  15]  were  available 
for  species  [1,  5]  of  other  nonicosahedral  symmetries, 
and  these  produced  radius-capacitance  points  that  were 
somewhat  off  the  regression  lines  shown  in  the  inset  of 
Fig.  1  and  in  Fig.  3.  However,  the  available  data  for 
points  corresponding  to  species  of  these  other  symme¬ 
tries  were  too  sparse  to  determine  scaling  lines  for  them. 
Thus,  we  only  display  here  the  data  and  points  for  C2 
and  V-2  nonicosahedral  fullerenes. 

The  DFT  I  value,  7.70  eV,  that  is  calculated  for  Ceo 
and  presented  in  part  A  of  Table  I  is  in  fairly  good  agree¬ 
ment  with  the  experimental  value  [4]  of  7.57  eV.  How¬ 
ever,  the  DFT  A  value  of  2.90  eV  for  Cgo  differs  more 
significantly  from  the  experimental  value  [3]  of  2.65  eV. 
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FIG.  2.  Icosahedral  fullerene  ionization  potentials  (black  cir¬ 
cles)  and  electron  affinities  (black  squares),  as  calculated  in 
this  work  using  density  functional  theory,  are  plotted  versus 
reciprocals  of  the  n-carbon  neutral  molecules’  average  radii, 
1  /rn.  Values  are  from  part  A  of  Table  I.  The  solid,  curved  line 
for  ionization  potential  function  I  and  dashed,  curved  line  for 
electron  affinity  function  A  are  determined  via  second-order 
regression,  for  which  equations  and  parameters  appear  in  the 
figure.  Dotted  straight  lines  present  scaling  of  /  and  A  ac¬ 
cording  to  a  purely  classical  electrostatic  model.  Points  rep¬ 
resenting  the  carbon  atom  would  lie  on  the  solid  and  dashed 
lines  for  I  and  A,  were  the  axes  extended.  See  text. 


This  error  is  believed  to  arise  primarily  from  the  self¬ 
interaction  in  the  DFT  local  density  approximation.  This 
self-interaction  might  be  removed,  but  only  with  con¬ 
siderable  computational  cost.  Using  a  different  density 
functional,  the  effect  of  the  self-interaction  was  studied 
for  atoms  and  small,  well-behaved  molecules  [21].  From 
these  considerations,  the  average  effect  of  not  correct¬ 
ing  I  for  the  self-interaction  is  expected  to  be  negligible, 
while  it  would  be  expected  to  make  A  values  from  DFT 
too  large  by  as  much  as  0.2  to  0.3  eV,  consistent  with 
what  is  observed  here  for  Ceo- 

In  Fig.  1,  the  solid  line  and  the  dashed  line  in  the  large 
graph  are  fit  by  linear  regression  to  the  icosahedral  and 
nonicosahedral  ( rn,Cn )  points  from  parts  A  and  B,  re¬ 
spectively,  of  Table  I.  The  inset  within  the  figure  provides 
a  more  detailed  view  of  the  nonicosahedral  capacitance 
scaling  line  and  the  points  that  define  it,  all  of  which  rep¬ 
resent  fullerenes  of  C2  or  V 2  symmetry  types,  as  indicated 
in  the  table.  Each  ( fn,Cn )  point  on  the  nonicosahedral 
scaling  line  represents  values  associated  with  the  lowest 
energy  [19,  20]  neutral  geometry  for  that  value  of  n. 

In  the  inset  of  Fig.  1,  it  also  is  seen  that  two  positions 
are  possible  for  the  radius-capacitance  point  of  C76.  This 
is  because  two  experimental  values  are  available  for  Re- 
Boltalina  et  al.  [4],  the  source  of  the  In  values  for  the 
other  nonicosahedral  fullerenes,  measured  Itq  =  7.34  eV, 


FIG.  3.  Nonicosahedral  fullerene  ionization  potentials  (black 
circles)  and  electron  affinities  (black  squares)  are  plotted  ver¬ 
sus  reciprocals  of  the  n-carbon  neutral  molecules’  average 
radii,  l/rn.  Values  of  In  and  An  are  from  experimental  re¬ 
sults  reported  in  the  literature  and  presented  in  part  B  of 
Table  I.  Values  of  l/rn  are  from  theory  [14,  15].  As  in  Fig.  2, 
the  solid,  curved  line  for  ionization  potential  function  I  and 
dashed,  curved  line  for  electron  affinity  function  A  are  deter¬ 
mined  via  second-order  regression,  for  which  equations  and 
parameters  appear  in  the  figure.  Unlike  the  nonlinear  regres¬ 
sion  lines  in  Fig.  2,  above,  the  lines’  intercepts  at  (l/rn)  =  0 
were  fixed  in  advance  at  the  value  of  the  graphene  work  func¬ 
tion,  4.7  eV.  See  text. 


while  Steger  et  al.  [2]  measured  176  =  7.10  eV.  Calculating 
C76  with  the  former  value  leads  to  a  radius-capacitance 
point  (black  square)  that  lies  off  the  capacitance  scaling 
line  defined  by  points  for  the  other  Ci  species  and  2?2 
species.  On  the  other  hand,  using  the  Steger  et  al.  value 
for  It 6  produces  a  point  (open  square)  that  is  closer  to 
the  line.  Thus,  we  judge  the  Steger  /76  value  and  the 
radius-capacitance  point  it  produces  as  more  consistent 
and  probably  more  accurate  than  that  of  Boltalina.  How¬ 
ever,  due  to  the  discrepancy,  we  did  not  use  points  for 
C76  in  any  of  the  regressions  in  this  paper. 

The  linear  fits  that  define  the  icosahedral  and  nonicosa¬ 
hedral  capacitance  scaling  lines  both  are  very  strong,  as 
seen  from  the  large  values  of  R 2  for  each  displayed  in 
Fig.  1.  Additionally,  it  is  seen  there  that  the  icosahedral 
and  nonicosahedral  scaling  lines  nearly  intersect  at  the 
Ceo  point.  This  might  be  explained  by  observing  that 
Ceo  is  a  kind  of  “progenitor”  structure  for  both  symme¬ 
try  types,  because  it  is  the  only  fullerene  whose  points 
all  rest  on  a  single  sphere.  As  more  carbon  atoms  are 
added  to  produce  larger  fullerenes,  two  growth  paths  are 
possible:  either  the  carbons  are  added  in  such  a  way 
as  to  maintain  the  high-symmetry,  truncated  icosahedral 
structure  of  C6o,  or  they  are  added  in  such  a  way  as  to 
reduce  the  symmetry  of  the  molecule. 
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FIG.  4.  Sketches  for  three  basic  cases,  (a)-(c),  illustrating 
the  relationship  between  the  scaling  versus  l/rn  of  the  elec¬ 
tron  detachment  energies  /  and  A  (depicted  at  left)  and  the 
scaling  versus  r„  of  the  quantum  capacitances  Cn  derived 
from  those  detachment  energies  (depicted  at  right).  Specif¬ 
ically,  the  sketches  show  the  manner  in  which  the  intrinsic 
nonlinearities  seen  in  the  quantum  electron  detachment  en¬ 
ergy  graphs  at  left  in  parts  (a)  and  (c)  of  the  figure,  as  well  as 
in  Figs.  2  and  3,  are  associated  with  the  nonclassical,  nonzero 
intercepts  shown  in  the  corresponding  capacitance  graphs  at 
right,  as  well  as  in  Fig.  1.  By  contrast,  as  shown  in  part 
(b)  above,  purely  classical  linear  scaling  of  the  detachment 
energies  leads  the  capacitance  intercept  to  vanish.  Classical 
scaling  is  not  realized  for  the  fullerenes,  but  is  a  very  good 
approximation  for  molecular  wires  [7].  See  text. 


We  observe  further  in  Fig.  1  that  the  radius- 
capacitance  point  for  the  C  atom  [6]  lies  right  on  the 
scaling  line  defined  by  the  ( rn,Cn )  points  for  the  icosa¬ 
hedral  fullerenes.  In  fact,  if  the  C  atom  point  were  to 
be  included  in  the  regression  with  the  points  for  the 
icosahedral  fullerenes,  the  correlation  coefficient  would 
increase  to  R2  =  0.9994.  For  this  reason,  the  carbon 
atom  may  be  thought  of  as  a  progenitor  capacitive  struc¬ 
ture  for  the  icosahedral  fullerene  capacitors;  the  icosahe¬ 
dral  molecules  share  a  certain  similarity  to  the  atom  in 
their  valence  energetics. 

The  regression  lines  in  Fig.  1  have  equations  of  the 


TABLE  I.  Quantum  capacitances  Cn  for  n-carbon  fullerenes 
tabulated  as  a  function  of  their  average  radii  rn.  Values  of 
Cn  are  calculated  via  Eq.  (1)  from  each  fullerene’s  ionization 
potential  In  and  electron  affinity  A„.  In  part  A  of  the  table,  In 
and  An  are  determined  for  icosahedral  fullerenes  via  detailed 
density  functional  theory  calculations  performed  in  this  work, 
while  in  part  B  they  are  from  experiments  on  nonicosahedral 
fullerenes  by  Boltalina  et  al.  [3,  4],  except  as  noted.  Data  are 
tabulated  here  only  for  nonicosahedral  fullerenes  of  C2  and  T> 2 
symmetries.  See  text. 


Quantum 

No.  of 

Ioniz- 

Capacit- 

Carbon 

Avg. 

ation 

Electron 

ance  fm. 

Sym- 

Atoms, 

Radius®, 

Potl., 

Affinity, 

Eq.  (1)  b, 

metryc 

n 

rn 

In 

An 

Cn 

(A) 

(eV) 

(eV) 

(+e/V) 

A.  Icosahedral  Fullerenes 

60 

3.548 

7.6996 

2.9005 

0.2084 

Th 

180 

6.135 

6.9189 

3.4520 

0.2884 

Th 

240 

7.073 

6.5460 

3.5379 

0.3324 

Th 

320d 

8.038 

5.7332 

4.1975 

0.6512 

Th 

540 

10.553 

6.0483 

3.9212 

0.4701 

Th 

720 

12.166 

5.9342 

4.0278 

0.5245 

Th 

960 

14.034 

5.7617 

4.1097 

0.6053 

Th 

1500 

17.522 

5.5736 

4.2276 

0.7430 

Th 

2160 

21.014 

5.4419 

4.3068 

0.8810 

Th 

B. 

Nonicosahedral  Fullerenes 

76® 

3.991 

7.10 

2.89 

0.238 

V2 

80 

4.094 

7.30 

3.17 

0.242 

V2 

82 

4.142 

7.25 

3.14 

0.243 

c2 

84 

4.193 

7.17 

3.14 

0.248 

V2 

86 

4.241 

7.16 

3.23 

0.254 

c2 

90 

4.339 

7.09 

3.27 

0.262 

c2 

94 

4.432 

6.96 

3.21 

0.267 

£2 

96 

4.482 

6.92 

3.28 

0.275 

V  2 

98 

4.53 

6.95 

3.26 

0.271 

c2 

100 

4.572 

6.95 

3.32 

0.275 

V  2 

a  Average  radius  rn 

is  the  arithmetic  mean 

of  distances  of  all  the 

carbon  atoms  in  a 

fullerene  from  a  central  point  within  the 

molecule,  as  determined  from  icosahedral  fullerene  geometries 
optimized  in  this  work  and  from  nonicosahedral  geometries  due 
to  Yoshida  and  Osawa  [14,  15]. 

°  Quantum  capacitances  Cn  are  reported  in  +e/V,  fundamental 
units  of  positive  charge  per  Volt.  Multiplication  by 
1.602188  X  10— 19  Coulombs  per  fundamental  unit  of  charge 
converts  these  capacitances  to  the  more  familiar  mks  units  of 
Farads.  Using  these  capacitance  units,  the  permittivity  of  free 
space  is  e0  =  5.526350  x  10“3  +e/V-A. 
c  For  nonicosahedral  fullerenes,  symmetry  type  is  taken  from 
Shao  and  Gao  [19,  20]  for  lowest  energy  fullerene  with  that 
number  of  carbon  atoms. 

d  C320  has  an  open-shell  electronic  structure,  unlike  that  of  the 
other  icosahedral  fullerenes  listed,  which  are  closed  shell. 
e  C76  point  not  used  in  regressions.  Value  of  J-jq  from  Steger  et 
al.  [2].  See  text. 


form 


C  =  4:ne0nr  +  C0,  (2) 

where  e0  =  5.526350  x  10-3  +e/V— A  is  the  permittivity 
of  free  space.  The  first  term  on  the  right  of  Eq.  (2)  is 
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Faraday’s  and  Maxwell’s  law  for  spherical  capacitors  in 
classical  electrostatics  [22,  23].  Thus,  the  dimensionless 
parameter  k,  which  determines  the  slope  of  the  scaling 
line,  is  the  analog  of  a  classical  dielectric  constant,  but 
for  individual  fullerene  molecules.  Equation  (2)  implies 
the  values  0.562  and  1.043  for  k  in  the  icosahedral  and 
nonicosahedral  cases,  respectively.  A  strikingly  nonclas- 
sical  [6]  feature,  though,  is  the  additional  constant  term 
Cq  that  is  required  to  account  for  the  empirically  ob¬ 
served  nonzero  values  of  the  capacitance  intercepts  at 
r  =  0  . 

In  Fig.  2,  the  icosahedral  fullerenes’  electron  detach¬ 
ment  energies  from  part  A  of  the  table  are  plotted  ver¬ 
sus  the  reciprocal  of  the  icosahedral  fullerenes’  average 
radii.  Then,  two  second-order  fits  determine  the  regres¬ 
sion  curves  through  the  (l/rn,  /„)  and  (l/r„,  An)  points, 
respectively.  The  values  of  7„  and  An  are  seen  to  progress 
smoothly  along  these  curves,  from  those  for  Cgo  at  the  far 
right  to  very  nearly  equal  intercepts  7^  «  A^  «  4.7  eV  at 
l/r„  =  0  (the  limit  of  very  large,  macroscopic  fullerenes). 
There,  the  energy  gap  ( In  —  An)  vanishes  to  the  accuracy 
of  the  calculations.  Such  an  intercept  sometimes  is  taken 
to  be  an  estimator  for  the  work  function  of  graphene  [4] , 
which  the  carbon  surface  of  a  very  large  fullerene  locally 
resembles.  This  is  because,  in  the  limit  of  very  large 
n  and  rn,  a  fullerene’s  surface  is  effectively  flat,  locally, 
and  consists  of  almost  purely  hexagonal  arrangements  of 
atoms,  since  it  still  contains  only  a  small,  fixed  number  of 
pentagonal  “impurities”  (12  for  any  fullerene).  The  limit¬ 
ing  value  of  4.7  eV  determined  here  is  in  close  agreement 
with  prior  estimates  of  4.6  eV  for  the  work  function  of 
graphene  [24]  and  graphite  [25] ,  which  are  thought  to  be 
approximately  equal. 

Fig.  3  analogously  plots  data  from  part  B  of  Table  I 
for  the  nonicosahedral  fullerenes.  A  difference  in  the  fit¬ 
ting  procedure  for  Fig.  3,  though,  is  that  we  fixed  the 
intercepts  with  the  energy  axis  at  the  value  of  the  work 
function  of  graphene:  7^  =  Aoc  =  4.7  eV.  Then,  we  per¬ 
formed  the  regression  analyses  to  determine  the  curves. 

The  curves  in  Figs.  2  and  3  are  fit  with  high  confi¬ 
dence  to  the  points  via  second-order  regression,  defining 
polynomial  expansions: 

In  =  loo  +  Xi(l/rn)  +  2TliXlrn)2  (3a) 

An  =  Aqq  +  A^(l/r„)  +  -Ty^l/fn)2.  (3b) 

Parameters  Xk  and  Tfc  represent  the  slopes  and  cur¬ 
vatures,  respectively,  of  the  detachment  energy  scaling 
curves.  The  fact  that  the  progression  of  the  detachment 
energies  from  those  of  Ceo  to  the  limiting  values,  7 <*,  and 
Aoo,  occurs  along  second-order,  curved  functions  of  1/r 
contrasts  with  the  linear  functions  of  1/r  anticipated  by 
simple,  classical  electrostatic  models  [10,  11],  as  noted 
above. 

Such  linear  functions  are  displayed  in  Fig.  2  as  the  dot¬ 
ted  lines  that  were  constructed  to  connect  the  76 0  point 
and  the  A60  point  to  7^  and  A <*,,  respectively.  It  is  seen 


that  these  classical  lines  give  reasonable  approximations 
for  all  the  I  and  A  values  determined  by  DFT  for  the 
icosahedral  fullerenes.  A  similarly  constructed  classical 
linear  approximation  would  be  even  closer  to  the  solid 
and  dashed  scaling  curves  shown  in  Fig.  3  for  the  non¬ 
icosahedral  fullerene  detachment  energies.  In  fact,  the 
dotted  lines  are  not  drawn  in  for  the  classical  approx¬ 
imation  because  they  would  be  too  close  to  the  curves 
determined  by  the  experimental  values,  and  would  ob¬ 
scure  them. 

The  classically  expected  lines  in  Fig.  2  have  the  equa¬ 
tions  7=10. 642(l/r) +4.7  and  A=— 6.385(l/r)+4.7.  Val¬ 
ues  of  In  and  An  from  along  these  classically  expected 
detachment  energy  functions  are  used  to  determine  the 
classically  expected  capacitances  for  each  f  =  f„,  and 
thereby  plot  the  points  (open  circles)  and  their  dotted 
classical  capacitance  regression  line  in  Fig.  1,  which  has 
an  intercept  Co  that  vanishes  to  the  accuracy  of  the  data. 

By  definition,  the  dotted  classical  capacitance  regres¬ 
sion  line  in  Fig.  1  goes  through  the  radius-capacitance 
point  for  Ceo-  Even  so,  it  is  seen  that  it  does  not  give 
good  estimates  for  the  capacitances  of  the  other  icosa¬ 
hedral  fullerenes  (along  the  solid  regression  line).  How¬ 
ever,  the  dotted  classical  line  provides  good  quantita¬ 
tive  estimates  for  the  capacitances  of  the  nonicosahe¬ 
dral  fullerenes  (along  the  dashed  regression  line).  This 
is  made  particularly  clear  by  the  close  proximity  of  the 
two  scaling  lines  in  the  inset  of  the  figure. 

In  order  to  preserve  the  detail  presently  seen  in  Fig.  2, 
the  points  for  a  single  carbon  atom  are  not  included  on 
the  graph  there.  Nonetheless,  just  as  the  point  repre¬ 
senting  the  carbon  atom  or  Ci  lies  on  the  quantum  ca¬ 
pacitance  scaling  line  for  the  icosahedral  fullerenes  in 
Fig.  1,  the  points  representing  Ci  (with  l/rq  =  1.087 
A-1,  7=11.26  eV,  and  A  =  1.26  eV)  would  be  seen  to  lie 
on  the  icosahedral  fullerene  quantum  scaling  curves  for  7 
and  A  in  Fig.  2,  were  the  scales  of  the  axes  in  the  figure 
extended.  Inclusion  of  these  Ci  points  in  the  regressions 
does  induce  a  small  energy  gap  of  approximately  0.3  eV 
between  7  and  A  in  the  plot  at  1/r  =  0.  However,  fits 
versus  1/r  of  fullerene  detachment  energies,  including  the 
Ci  points,  both  are  very  strong,  with  R2  =0.997  for  both 
7  and  A.  Further,  if  the  induced  gap  is  eliminated  from 
these  fits,  by  constraining  both  intercepts  at  1/r  =  0  to 
assume  the  value  of  the  graphene  work  function,  4.7  eV, 
the  fits  continue  to  be  very  strong,  with  R 2  >  0.992  for 
both  curves.  Thus,  an  additional  important  result  of  this 
work  is  that  both  the  quantum  capacitances  and  the  elec¬ 
tron  detachment  energies  of  icosahedral  fullerenes  scale 
smoothly  from  the  limit  of  a  single  carbon  atom,  through 
points  for  Ceo,  to  the  points  representing  a  graphene-like 
structure,  in  the  limit  of  n  =  oo  at  1/r  =  0. 


III.  DISCUSSION  AND  ANALYSIS 

Analysis  of  the  foregoing  results  reveals  that  the 
nonzero  quantum  capacitance  intercepts  Co  for  the 
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fullerenes,  as  seen  in  Fig.  1,  are  a  consequence  of  the  non- 
classical,  nonlinear  scaling  of  the  electron  detachment  en¬ 
ergy  functions  in  Figs.  2  and  3,  and  vice  versa.  To  demon¬ 
strate  this  algebraically,  we  substitute  Eqs.  (3)  in  Eq.  (1), 
then  expand  the  result  in  a  geometric  series  through  first 
order  in  rn: 


CL  = 


(  1  > 


1 


TA  -  T/ 


2  L (A/  —  Aa)2J 


(4) 


To  derive  this  result,  we  assume  that  ~  Aoo,  so 
that  their  difference  vanishes,  at  least  approximately,  as 
discussed  above.  Then,  order-by-orcler  comparison  of 
Eq.  (4)  to  Eq.  (2)  enables  us  to  evaluate  in  terms  of 
the  slope  and  curvature  of  the  detachment  energy  scal¬ 
ing  curves  the  molecular  dielectric  constant 


1 


K  = 


(  1  ) 

4"7r 6q  V  A/  —  Ayt  / 


and  the  capacitance  intercept 


1  r 


Cn=  - 


T A  -  Tj 


2  L  (A/  —  XA)2- 

=  87r2£g k2(ta  -  tj). 


(5) 

(6a) 

(6b) 


Equations  (4)  through  (6)  involve  no  classical  approxi¬ 
mations. 

Now,  applying  Eq.  (6b)  to  Fig.  2,  for  example,  it  is 
seen  that  the  curvatures,  ta  and  tj,  of  the  second-order 
regression  curves  for  In  and  An  are  of  opposite  sign,  so 
the  two  curvatures  do  not  cancel  in  Eq.  (6b)  and  Co  does 
not  vanish.  Thus,  from  Eq.  (6b),  the  positive  nonzero  ca¬ 
pacitance  intercept  may  be  regarded  as  a  consequence  of 
the  nonzero  curvatures  of  the  detachment  energy  scaling 
curves  plotted  in  Fig.  2.  This  case  also  is  summarized  in 
Fig.  4(a). 

Evaluating  Eq.  (6b)  for  the  case  depicted  in  Fig.  3 
once  again  yields  a  nonzero  value  of  Co-  Here,  however, 
the  difference  between  the  curvatures  is  negative,  because 
r/  >  0  and  ta  ~  0.  Thus,  C0  <  0,  as  plotted  in  Fig.  1 
and  depicted  schematically  in  Fig.  4(c). 

Conversely,  we  also  may  regard  the  nonzero  curva¬ 
tures  for  the  detachment  energy  graphs  as  a  consequence 
of  a  nonzero  capacitance  intercept.  Considering  only 
the  valence  ionization  potential  and  starting  from  the 
purely  classical  electrostatic  formulas  of  Smith  [10]  and 
of  Leach  [11],  we  may  write: 


In  —  loo  I/^Cn  (7a) 

<7b) 

The  second  line  in  the  previous  equation  follows  from 
substitution  of  quantum  capacitance  equation  (2)  into 
Eq.  (7a),  then  applying  a  first-order  geometric  series  ex¬ 
pansion  to  the  result.  Analogously,  using  Smith’s  [10] 
classical  expression  for  the  electron  affinity  as  a  starting 
point,  we  may  write: 


A  —  A 
~  Aoo 


1/2  Cn  (8a) 

(sb)/;+2C“(8i)2(/;)2-(8b) 


In  Eqs.  (7)  and  in  Eqs.  (8),  to  obtain  the  third  term 
on  the  right  that  corresponds  to  a  nonzero  curvature  like 
that  observed  in  the  quantum  results  plotted  in  Fig.  2,  it 
is  essential  that  the  expansions  above  utilize  both  terms 
of  Eq.  (2).  This  differs  from  the  purely  classical  expan¬ 
sion  procedure  employed  by  Smith  [10]  and  by  Leach  [11] 
that  utilizes  only  the  first  term  on  the  right  of  Eq.  (2), 
effectively  setting  Co  =  0  in  Eqs.  (7b)  and  (8b). 

The  expressions  in  Eq.  (7b)  and  Eq.  (8b)  for  the  coef¬ 
ficients  of  the  powers  of  1  jrn  can  be  verified  to  be  rea¬ 
sonably  accurate,  quantitatively.  For  example,  we  ob¬ 
tain  from  the  regression  equations  in  Fig.  1  the  values 
47re0«:  =  4.69  +e/V-A  and  Co  =  0.0582  +e/V,  in  the  case 
of  the  icosahedral  fullerenes.  Using  these  in  Eq.  (7b),  the 
values  for  the  first  and  second-order  coefficients  are  12.82 
eV-A  and  -19.13  eV-A2,  respectively.  These  compare  fa¬ 
vorably  to  the  respective  values  16.61  eV-A  and  -20.95 
eV-A2  seen  in  the  regression  equation  for  the  accurate 
ionization  potentials  in  Fig.  2. 

Continuing  the  derivation,  from  a  comparison  of 
Eqs.  (7b)  and  (8b)  with  Eqs.  (3a)  and  (3b),  respectively, 
one  may  obtain  approximate  relations  for  the  slopes  and 
curvatures  of  /  and  A  with  respect  to  1  /r: 


A  i  =  —  Aa  = 


8ireon 


(9) 


Ti  =  - TA  =  — 4C0A|.  (10) 

Equations  (7b)  through  (10)  incorporate  quantum  effects 
in  the  nonzero  curvatures,  but  can  only  be  approximate 
because  of  the  use  of  the  classical  expressions  in  Eqs.  (7a) 
and  (8a)  as  starting  points  in  their  derivations.  Nonethe¬ 
less,  Eqs.  (9)  and  (10),  along  with  Eq.  (3),  do  provide  a 
rationale  for  the  near  symmetry  of  I  and  A  as  a  function 
of  l/r  that  is  manifest  in  both  Fig.  2  and  Fig.  3. 

From  the  accurate  quantum  results  embodied  in  the 
regression  equations  within  these  figures  one  can  verify, 
as  well,  that  tj  «  —  ta-  In  the  values  of  the  regression 
parameters  for  the  icosahedral  fullerenes  given  in  Fig.  2, 
there  is  order  of  magnitude  agreement  of  77  and  ta  with 
Eq.  (10).  In  Fig.  3,  for  the  nonicosahedral  fullerenes, 
the  approximate  agreement  with  this  curvature  relation 
is  manifest  if  one  recognizes  there  that  ta  ~  0. 

Equation  (10)  also  shows  explicitly  that  a  nonzero  ca¬ 
pacitance  intercept  results  in  a  nonzero  curvature  for  I 
and  for  A  as  a  function  of  1  /rn,  as  asserted  above.  Still 
further,  from  that  equation  it  is  clear  that  the  sign  of  the 
capacitance  intercept  Co  and  the  signs  of  the  curvatures 
t  of  the  detachment  energy  scaling  graphs  are  mutually 
related.  Thus,  for  example,  just  as  a  positive  value  of  Co 
for  the  icosahedral  fullerenes  is  ensured  by  (and  ensures) 
a  negative  value  of  77,  a  negative  value  of  Co  in  the  case 
of  the  nonicosahedral  fullerenes  is  associated  with  a  pos¬ 
itive  value  of  tj.  This  can  be  verified  empirically  from 
the  results  in  Figs.  1,  2,  and  3. 

More  generally,  in  Fig.  4,  we  map  out  schematically 
three  major  cases  of  this  relationship  between  the  quan- 
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turn  capacitance  intercepts  and  the  curvatures  in  de¬ 
tachment  energy  scaling  graphs.  These  apply  for  the 
fullerenes  studied  here  and  for  molecular  wires  studied  in 
prior  work  [7].  Other  cases  may  arise  for  still  other  quan¬ 
tum  systems.  For  the  cases  depicted  in  Figs.  2,  3,  and  4, 
the  near  symmetry  in  the  scaling  of  /  and  A  about  the 
work  function  in  the  fullerenes  strongly  resembles  that 
observed  previously  for  clusters  of  metal  atoms  [13,  26] 
and  affirms  their  classification  as  semimetals. 

Considering  the  results  displayed  in  Fig.  1  from  a  more 
qualitative  viewpoint,  the  two  different  scaling  lines  for 
the  differently  shaped  icosahedral  and  nonicosahedral 
fullerenes,  respectively,  are  consistent  with  the  classical 
notion  that  capacitances  depend  strongly  on  the  shapes 
of  capacitors.  In  this  connection,  we  observe  further  that 
the  nonicosahedral  fullerene  capacitances  lie  above  the 
scaling  line  for  the  icosahedral  fullerenes.  This  is  consis¬ 
tent  with  the  facts  that  (a)  isoperimetric  principles  [27] 
of  classical  electrostatics  indicate  that  the  capacitance  of 
a  conductor  of  any  shape  is  proportional  to  the  square 
root  of  its  surface  area,  and  that  (b)  the  more  aspherical 
nonicosahedral  species  should  have  larger  surface  areas 
for  a  given  average  radius  than  do  icosahedral  fullerenes. 

In  addition,  we  see  in  Fig.  1  that  the  capacitance-radius 
point  for  open-shell  C320  (the  x )  lies  well  above  the  scal¬ 
ing  line  for  the  other  icosahedral  fullerenes,  which  have 
closed-shell  valence  electron  configurations.  This  outlier 
point  suggests  that  electronic  structure  (esp.,  the  mul¬ 
tiplicity),  as  well  as  geometric  structure  of  the  nuclear 
framework,  plays  an  important  role  in  determining  the 
effective  shape  and  dimensions  of  a  fullerene  capacitor. 
A  related  explanation  of  the  inordinately  large  capaci¬ 
tance  of  C320  is  that  its  open-shell  electronic  structure 
is  more  diffuse.  This  produces  a  larger  effective  radius 
and  surface  area  of  its  valence  electron  distribution  for 
its  value  of  fn  than  would  a  closed-shell  electron  distri¬ 
bution. 

IV.  SUMMARY  AND  CONCLUSIONS 

In  summary,  we  have  studied  the  scaling  of  fullerene 
quantum  capacitances  as  a  function  of  their  average  radii 
fn  and  the  scaling  of  the  molecules’  valence  electron  de¬ 
tachment  energies  as  a  function  of  l/rn.  Detailed  density 
functional  theory  calculations  of  electron  detachment  en¬ 
ergies  for  icosahedral  fullerenes  and  results  from  earlier 
experimental  measurements  of  detachment  energies  for 
nonicosahedral  fullerenes  were  employed  for  this  purpose. 
Linear  scaling  was  found  for  the  capacitances  of  both 
symmetry  types;  however,  this  scaling  occurred  along  two 
different  scaling  lines,  one  for  icosahedral  fullerenes  and 
one  for  nonicosahedral  fullerenes  of  T> 2  and  C2  symme¬ 
tries.  (See  Fig.  1.) 

Though  linear  scaling  was  found  for  the  capacitances 
as  a  function  of  rn,  nonclassical,  no?rlinear  scaling  was 
found  for  the  fullerene  detachment  energies  as  a  function 
of  l/rra.  These  nonlinear  trends  were  seen  to  resemble 


those  for  clusters  of  metal  atoms.  Further,  a  proof  in  Sec¬ 
tion  III  shows  that  the  nonlinear  behavior  seen  in  Figs.  2 
and  3  for  the  detachment  energies  leads  to  the  nonclassi¬ 
cal,  nonzero  intercepts  seen  in  Fig.  1  for  the  capacitance 
scaling  lines,  and  conversely. 

Thus,  classical  scaling  is  not  realized  in  the  fullerenes. 
Classical  models  commonly  employed  to  represent  their 
detachment  energies  and  capacitances  yield  reasonably 
good  quantitative  estimates,  but  fail  to  reproduce  key 
features  of  the  accurate  quantum  scaling  trends.  Espe¬ 
cially,  they  miss  the  nonlinearity  of  the  detachment  ener¬ 
gies  as  a  function  of  1  /rn  and  nonzero  intercepts  for  the 
capacitance  as  a  function  of  fn.  They  also  fail  to  account 
for  the  differences  between  the  trends  for  the  icosahedral 
and  nonicosahedral  species.  Instead,  a  simple  set  of  alge¬ 
braic  equations  is  derived  in  Section  III  to  explain  these 
quantum  behaviors  and  the  relationships  among  them, 
as  well  as  to  contrast  them  with  classically  expected  be¬ 
haviors. 

Quantum  capacitance  scaling  of  the  icosahedral  and 
nonicosahedral  fullerenes  also  is  contrasted  with  classi¬ 
cally  expected  capacitance  scaling  in  Fig.  1,  while  the 
departures  from  classical  detachment  energy  scaling  are 
shown  in  Figs.  2  and  3.  The  relationships  discovered  be¬ 
tween  the  nonclassical  detachment  energy  scaling  behav¬ 
iors  and  the  nonclassical  features  in  quantum  capacitance 
scaling  are  depicted  conceptually  in  Fig.  4. 

Another  key  finding  of  this  work  is  that  points  repre¬ 
senting  the  carbon  atom  lie  on  scaling  lines  for  icosahe¬ 
dral  fullerenes,  as  do  points  for  C6o>  while  the  detachment 
energy  scaling  curves  also  intersect  points  representing  a 
graphene-like  structure  in  the  limit  of  very  large  radius. 
Thus,  the  icosahedral  fullerene  quantum  capacitances 
and  electron  detachment  energies  all  scale  smoothly  from 
a  single  carbon  atom  to  the  graphene  limit. 

Lastly,  the  algebraic  formalism  derived  to  show  the 
connection  between  nonclassical  behaviors  in  detach¬ 
ment  energy  scaling  and  quantum  capacitance  scaling 
for  fullerenes  should  be  useful,  as  well,  for  interpret¬ 
ing  quantum  scaling  trends  for  other  systems,  such  as 
atoms  [6]  and  diatomic  molecules  [8] .  Similarly,  the  equa¬ 
tions  and  conceptual  insights  developed  here  may  be  use¬ 
ful  in  providing  more  accurate,  simple  estimators  of  elec¬ 
tron  detachment  energies  for  other  homologous  series  of 
molecules  and  for  still  other  types  of  nanostructures. 
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